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We consider the collision of a dark soliton with an obstacle in a quasi-one- dimensional Bose 
condensate. We show that in many respects the soliton behaves as an effective classical particle of 
mass twice the mass of a bare particle, evolving in an effective potential which is a convolution of 
the actual potential describing the obstacle. Radiative effects beyond this approximation are also 
taken into account. The emitted waves are shown to form two counterpropagating wave packets, 
both moving at the speed of sound. We determine, at leading order, the total amount of radiation 
emitted during the collision and compute the acceleration of the soliton due to the collisional 
^ . process. It is found that the radiative process is quenched when the velocity of the soliton reaches 
in ' the velocity of sound in the system. 
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1 Introduction 



One of the many interesting aspects of the physics of Bose-Einstein condensation of ultracold 
atomic vapors is to open opportunities of studying mesoscopiclike phenomena in new types of 
setups. The advances in the production and propagation of Bose-Einstein condensates in more and 
more elaborate waveguides (magnetic or optical, microfabricated or not [1]) opens up the prospect 
of studying a rich variety of quantum transport phenomena for these intrinsically phase-coherent, 
finite-sized systems. In particular it has been possible to study quantum interference effects [2], 
Bloch oscillations and Landau-Zener tunneling [3], Josephson junctions [4], and superfluidity [5]. 

Pushing further the analogy in transport properties of mesoscopic systems and Bose-condensed 
vapors, one notices that, whereas in mesoscopic physics interaction effects are often difficult to 
understand, in Bose-Einstein condensates they are more easily accessible to theoretical description 
and have the advantage of covering a wide range of regimes, ranging from almost noninteracting 
atom lasers to strongly correlated systems. Along this line, the existence of nonlinearity in the wave 
equation, resulting in the existence of bright [6] a dark [7] solitons, appears as a natural - and rather 
simply understood - consequence of interaction on transport phenomena of quasi-one-dimensional 
Bose-condensed systems. 

In the present work we address the problem of transport of a dark soliton in a quasi-one- 
dimensional Bose-Einstein condensate. More precisely, we consider a guided Bose-Einstein conden- 
sate and theoretically study the propagation of a dark soliton encountering an obstacle on its way. 
In the appropriate limit [see Eq. Q below] the system is described by a one-dimensional nonlinear 
Schrodinger equation. This equation admits bright and dark solitonic solutions, depending on the 
sign of the interparticle interaction. The obstacle is modeled via an external potential, and this 
could correspond to different physical realizations, such as a heavy impurity, a (red or blue) detuned 
laser beam crossing the atomic beam, a bend, a twist, or a constriction in the shape of the guide. 

A soliton under the influence of a perturbation (here, the obstacle) sees its shape and velocity 
modified and may also radiate energy (see, e.g., Ref. [8]). Despite their mutual dependence, these 
two phenomena are not easily treated on the same theoretical footing. The evolution of the soliton's 
parameters is typically studied within the adiabatic approximation (see Ref. [9] and references 
therein), whereas radiative effects are not so easily described, because their influence on the soliton's 
parameters only appears at second order in perturbation theory (see the discussion in Section f4.4|) . 
However, it has been possible to treat both phenomena concomitantly in the case of bright solitons 
[8-12]. Concerning dark solitons, several studies of adiabatic dynamics have appeared [13-19], but 
until recently radiative effects have been treated mainly numerically [20-22]. 

In the present paper we study the dynamics of a dark soliton via perturbation theory. This 
method, based on the theory of linear partial differential equations, has been established in the case 
of the nonlinear Schrodinger equation with repulsive interaction in Refs. [23,24] (see also the earlier 
attempt [25]). Although our first interest lies in the physics of guided Bose-Einstein condensates, 
the method employed and the results displayed also apply to optical waveguides described by a 
one-dimensional (ID) nonlinear defocussing Schrodinger equation. 
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The paper is organized as follows. In Section [21 we present the basic ingredients of the model 
and the resulting equation governing the time evolution of the condensate wave function. In the 
framework of perturbation theory we then derive the equations determining the dynamics of the 
soliton and of the radiated part (Section |3J) . The results are analysed in Section |1J We show that 
one can devise a quite successful approximation that we denote as "effective potential theory," 
where the soliton is assimilated to a classical particle of mass twice the mass of a bare particle, 
evolving in an effective potential (Section I4.1|) . The agreement of this approximation with the 
results of the adiabatic approximation is verified even for a fine quantity such as the position 
shift induced on the trajectory of the soliton by the obstacle (Section 14, 2j) . We then consider in 
Section ^, 3l the radiated part and show that it is formed of backward- and forward-emitted phonons, 
which form two counterpropagating wavepackets mooving at the speed of sound. In the limit of 
large soliton's velocity we furthermore obtain in Section [4.41 an analytical expression for the total 
amount of radiation emitted by the soliton during the collision. In addition we show that (within our 
leading-order evaluation) a soliton reaching the velocity of sound does not radiate, and we propose 
a physical interpretation for this phenomenon. Finally we present our conclusions in Section 03 
Some technical points are given in the Appendixes. In Appendix A we recall the main properties 
of the spectrum of the operator governing the wave dynamics of the system around the solitonic 
solution. In Appendix B we briefly present the Lagrangian approach for deriving the dynamics of 
the parameters of a dark soliton. In Appendix C we show how to compute some integrals involved 
in the evaluation of the total amount of radiation emitted by the soliton. 



2 The model 

We consider a condensate confined in a guide of axis z and denote by n(z, t) the ID density of the 
system. The condensate is formed by atoms of mass m which interact via a two-body potential 
characterized by its 3D s-wave scattering length a sc . We consider the case of a repulsive effective 
interaction-i.e., a sc > 0. The condensate is confined in the transverse direction by an harmonic 
potential of pulsation u±. The transverse confinement is characterized by the harmonic oscillator 
length a± = (l/mwi) 1 ' 2 . 

With n 1D denoting a typical order of magnitude of n{z^ £), we restrict ourselves to a density 
range such that 

(a sc /a±) 2 < n w a sc < 1 • (1) 

This regime has been called "ID mean field" in Ref. [26]. In this range the wave function of 
the condensate can be factorized in a transverse and longitudinal part [27-29]. The transverse 
wave function is Gaussian (this is ensured by the condition n 1D a sc <C 1) and the longitudinal one, 
denoted by ip(z,t) [such that n(z, t) = \tp(z,t)\ 2 ], satisfies an effective ID Gross-Pitaevskii equation 
(see, e.g., [27-29]): 

Ti 2 ( 1 

—z— i>zz + \ U(z) + 2hiu ± a sc #| 2 \i/) = ihi/H. (2) 
1ml ) 

In Eq. ©, U(z) represents the effect of the obstacle. We restrict ourselves to the case of localized 
obstacle such that linxz-^oo U(z) = 0. Hence, we can consider that the stationary solutions of 



3 



Eq.© have at infinity an asymptotic density unperturbed by the obstacle. Besides, considering 
solutions without current at infinity, we impose the following form to the the stationary solutions: 

^sta(<M) = f(z) exp[-ifj,t/h] , with lim f(z) = y/n^ , (3) 

where rioo is the ID density far from the obstacle and fi = 2?iuj±_ a sc noo the chemical potential [30] . 

We note here that in Eq.© we have discarded very low densities in order to prevent the sys- 
tem from getting in the Tonks-Girardeau regime where the mean-field picture implicit in Eq. © 
breaks down [28, 31]. This can be intuitively understood as follows: it is natural to assume 
that the Gross-Pitaevskii scheme is valid-i.e., that the system can be described by a collective 
order parameter ^-only if the interparticle distance (of order n^ 1 ) is much smaller than the 
minimum distance £ over which ip can significantly vary [£ is the healing length, defined by 

Uoo) ]• The condition <C £ then imposes us to consider the 
regime a sc 3> (a sc /a±) 2 to which, from Eq. ((H). we restrict our study. If one considers, for in- 
stance, 87 Rb or 23 Na atoms in a guide with a transverse confinement characterized by wj_ = 2ir x 500 
Hz, the ratio a sc /a± is roughly of order 10 -2 and the restriction © still allows the density to vary 
over four orders of magnitude. 

In all the following we use dimensionless quantities: the energies are expressed in units of \i, 

— 1/2 

the lengths in units of £, and the time in units of h/fjL. ifi is also rescaled by a factor noo ; this 
corresponds to expressing the linear density in units of the density at infinity, n^. We keep the 
same notation z, t, U(z), and ip(z,t) for the rescaled quantities. Equation © now reads 

-\^ zz + {u{z) + \^\ 2 )i> = i^ t . (4) 

From Eq.©, the stationary solutions of Eq.Q are of type f{z) exp[— it], f being real, and a 
solution of 

~f zz + {u(z) + f 2 -l}f = 0, (5) 

with the asymptotic condition lim z _ + ± 00 /(;z) = 1. 

The method we will expose is quite general and applies to a broad range of potentials U (z), but 
for concreteness we will often display the explicit solutions of the problem in the case of a pointlike 
obstacle, where U(z) = X5(z); A > (< 0) corresponds to a repulsive (attractive) obstacle. For 
such an obstacle, the solution of Eq.© is 

/(*) = I ta ^n1 + 1 > n with a = - sinh^ ( *) . (6) 

Jy J \coth(|z|+a) if A<0 2 Vl A l/ 

In section 4 we will concentrate on perturbative aspects of the problem and consider the case 
of a weak potential U{z). For a pointlike obstacle, this corresponds to the limit |A| <C 1. In this 
case sinh _1 (2/|A|) ~ ln(4/|A|) and Eqs. © simplify to 

f(z)~l-~exp{-2\z\}. (7) 
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In the general case, one can design a simple treatment [29, 32] valid for any weak potential U(z) 
leading after linearization of Eq. © to the perturbative result f{z) = 1 + Sf(z) with 

i r+oo 

Sf(z)^~- dyU(y)exp{-2\z-y\} , (8) 

^ J —oo 

of which Eq.0 is a particular case. In section 4 it will reveal convenient to rewrite Eq.JHJ) in 
an other way: denoting by U(q) = f^dzU(z) exp(-iqz) the Fourier transform of U(z), one may 
equivalently express 5f defined in Eq.® as 

The stationary solutions of the problem being defined, let us now turn to the main subject of the 
present work and consider the case of time-dependent solutions corresponding to a dark soliton 
propagating in the system. The soliton will appear as a distortion of the stationary background, 
and it is here very natural to follow the approach of Frantzeskakis et al. [17] who write the wave 
function of the system as a product: 

V>(M) = 0(2, *)/(*) exp(-it) . (10) 

(j>(z,t) in Eq. ([TO]) accounts for the deformation of the stationary background f(z) exp(— it) caused 
by the motion of a soliton in the system. From Eq. (j3J we see that the unknown field <f)(z, t) is a 
solution of the following equation: 

i^ + \^z-{\<P\ 2 -l}^ = R[^} , (11) 

where 

«[0] = -^y + (/ 2 -l)(M 2 -l)0. (12) 

Far from the obstacle, f(z) = 1 and thus R{4>} = 0. In this case, the motion of a dark soliton 
in the system is described by the usual solitonic solution of the defocussing nonlinear Schrodinger 
equation [33] 

0(z,t) = <f>(z-Vt-b,9) , (13) 

where 

$(x, 9) = x( x i Q) + i V , with x(x, 0) = cos tanh(x cos 0) and sin# = y. (14) 

Equations (|13[) and (|14|) describe a dark soliton consisting in a density trough located at position 
Vt + b at time t. The phase change across the soliton is 29 — it. The choice of the parameter 
9 in [0, 7r/2] corresponds to a soliton moving from left to right with a velocity V = sin9 £ [0,1]. 
Note that a dark soliton has a velocity always lower than unity (which, in our rescaled units, is 
the velocity of sound [34]). When 9 = 0, the soliton is standing and its minimum density is zero; 
it is referred to as a black soliton. When 9^0 one speaks of a gray soliton. We display in Fig. 
1 the density profile and the phase of the wave function vp(z,t) [see Eq. pU[>] describing a soliton 
incident with velocity V = 0.4 on a repulsive point-like obstacle characterized by A = 0.5. 
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Figure 1: Upper plot: density profile of a 
dark soliton incident with velocity V = 0.4 
on a point-like repulsive obstacle U(z) = 
X8(z) (with A = 0.5). The arrow repre- 
sents the direction of propagation of the 
soliton. Lower plot: phase of the wave 
function ip(z, t = 0) describing the system. 
Across the soliton the phase of the wave 
function changes from tt — 6 to 8 (with 
V = sm6). 



3 Perturbation theory 

In the following we will set up the basis for a systematic perturbative expansion, and for properly 
identifying the orders of perturbation at which the expansion is done, it is customary to introduce 
an artificial multiplicative parameter e in the potential of the obstacle (otherwise of arbitrary form). 
We will see in the present section (and justify on physical grounds in the next one) that for an 
obstacle characterized by eU(z), the condition of small perturbation reads V 2 3> e U. Since the 
soliton velocity is always lower than unity (which is the speed of sound in our dimensionless units), 
this condition implies e(J « 1; i.e., Eqs. I)8|9|) hold. 

At initial times the soliton is unperturbed and described as in the previous section by (f)(z, t) = 
<&(z — Vt,8o) [$ is defined in Eq. IJ14JI ]: i.e., one considers a soliton incident from left infinity with 
velocity V = sin#o- The more important effect of the obstacle on the soliton is a modification of 
its shape; i.e., the parameters characterizing the soliton will become time dependent in the vicinity 
of the obstacle. Perturbation at next order describe the emission of radiations. One thus looks for 
solutions of Ea. (|llj) of the form 

<j)(z, t) = 4> BOl (z, z(t),9(t)) + 6<f>(z, t) , (15) 

where 

<f> ml (z, z(t),9(t)) = *(z - z(t),6(t)) , (16) 

describes a soliton which is characterized by the two parameters z(t) (decribing the center of the 
soliton) and 6(t) (describing the phase shift accross the soliton). 5(f) describes additional radiative 
components: 

5(j)(z,t) = e<M-M) + e 2 <MM) + - • (17) 

Equations (|15 |) -(|17 |) form the grounds of a secular perturbation theory where the time dependence 
of the parameters of the soliton permits to avoid the grow of secular perturbation in 5(f) (see, e.g., 
the discussion in Ref. [8]). 
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It is more appropriate to define 5(f) in Ea, ()15|) and the 0j's in Ea.p7|) as functions of z — z(t) 
than as functions of z. To this end, we define x = z — z(t) and choose to work with x and t as 
independent parameters rather than z and t. This corresponds to the transformation 

d z , d t -> d x , d t - id x . (18) 

Furthermore, in order to take into account the slow time dependence of the parameters of the soliton, 
it is customary to introduce multiple time scales: t n = e n t (n £ N). A time-dependent function 
could, for instance, depend on t via t\, indicating a weak time dependence (a ti dependence being 
related to an even weaker time dependence and a to dependence to a "normal" time dependence). 
Generically, time-dependent quantities will be considered as functions of all the t n 's, with 

d t = d t0 + ed tl +e 2 d t2 + ... . (19) 

In the following we will make an expansion at order e and it will suffice to consider only the fast 
time to and the first slow time t\. The soliton's parameters 9 and z are considered as functions 
0(ti) and z(t ,tt) [35]. 

Putting everything together, we see that, at order e, Eqs. (|15 )) -(|17 j) read explicitly 

<j)(z,t) = $(x,6(ti)) + e<f) 1 (x,t ,t 1 ) , with x = z - z(t Q , t x ) . (20) 

Equation(|ll|) is now rewritten taking the transformations (|18|) and (|19j) into account, with an 
expansion at order e. To this end, we have to take into account that R[4>] defined in Ea. ()12|) is a 
small quantity and can be written at first order in e as 

R[4>] ~ -d x $(x, 9) [d z 5f{z) + 2 9)5 f{z)] = e K(x, z) , (21) 

where z = x + z(tQ,t\) and 5f(z) is defined as in Eqs.© and ©, with an extra multiplicative 
factor e in U which has been written explicitly in the definition of 1Z on the right-hand-side (RHS) 
of Eq-d). 

We are now ready to expand Eq.tjTTj) in successive orders in e. The leading order reads 

-i^ + iz i0 $ :c + (|$| 2 -l)$ = 0, (22) 

implying that 

z to = sin 9 , (23) 

whence z can be written as 

z = t eane + zfa) , (24) 
where z{t\) is a still unknown function [36]. At next order in e one obtains 

idtofa = [-^d 2 x + is\n9d x + 2\^\ 2 -l\(f) 1 + ^ 2 (t)\+n-i9 tl ^e + izt 1 ^x ■ (25) 
Equation (|25|) can be rewritten as 

id t0 \(f)i) = H\(f)i) + a 3 \K) + iz tl \u e ) - -^-\n e ) , (26) 



7 



where = (0i,0*) T , \1Z) = (1Z,1Z*) T , 03 is the third Pauli matrice, and 

- ( -\d 2 x + ism9d x + 2\<$>\ 2 -I $ 2 \ . 

n ~\ -$* 2 £flg + isin00 a! -2|*| 2 + l J ' [Z<) 

7i is not diagonalizable, but can be put in a Jordan form in a manner similar to what has been 
done for the attractive nonlinear Schrodinger equation [37]. Its eigenfunctions and eigenvalues 
are presented in Appendix A In particular, \u e ) and \Q e ) appearing in Ea.H26[) belong to the 



generalized null space of TL; they verify H\io e ) = and 7i|O e ) = cos 2 #|cj e ). As well as its null 
space, 7i has two continuous branches of excitations which we denote by its "phonon spectrum." 
The corresponding eigenfunctions are denoted by |Sj^) with q £ R (see Appendix A I. 



It is physically intuitive that corresponding to the radiated part should be expanded over 
the phonon part of the spectrum of 7i: 

/+oo 
dqC^H)]^). (28) 
-oo 



A more technical argument for limiting the expansion (|28jl to the phonon components of the spec- 
trum of H is the following: one might think that a greater generality could be achieved by allowing 
|</>i) to have also components on \u> e ) and \£l e ), for instance. However, exactly as in the case of the 
bright soliton [12], these components can (and should) be imposed to remain zero for avoiding the 
appearance of secular terms in the evolution of the soliton's parameters. 



3.1 Evolution of the parameters of the soliton 

Applying (w e |o"3 and (O e |o"3 onto (|2T)|) and using the orthogonality relations (|A6|) . one obtains the 
equations of evolution of the parameters of the soliton: 



40, 



and 



tl cos 2 6 = -(u e \n) = -2Re ij dx$ x K*(x,x + z)\ 

= j Re dzR*[(j) sol } a 2 -0 sol | , (29) 



4z h cos 2 = — 5_— {Sl e \K) = -2 Re { f dx ^ e TZ*(x } x + z)\ 

l COS 9 [J-oo J 



-oo 



2 Re <j / dz R*[(t> sol ] d <p sol \ . (30) 



The set of equations l|23 |l .l|29 |) and ()3U|) describe the time evolution of the soliton's parameter. 
The same equations are obtained via adiabatic approximation which is a simpler variationnal ap- 
proximation where radiative effects are neglected [see Appendix B[ Eqs. ()B9|) and (|B1Q|) ] . This is 
evident in the case of Eq. (|29|) which is the slow time analogous to Eq. (|B9|) (since 9 = eO^). In a 
similar way, the prescription ((TO)) indicates that z = z to + ez tl ; combining Eqs. 1(2*3^) and (|30J) . one 
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sees that the equations of evolution of z obtained in the present section correspond to the multiple- 
time expansion of Eg. ()B10j) . As a side result of this exact correspondance of the time evolution of 
the soliton's parameters, we obtain here that, as in the adiabatic approach, sin[#(ti — ► ±00)] = V 
(see the discussion at the end of Appendix B ) and the quantity z appearing in Eq. ({21]) is identical 



when t\ — ► +00 to the one defined in Eq. (|48|) . 

A technical remark is in order here. One can notice that in Ea.(|15|) we did not consider the 
most general variational form for the solitonic component of the wave function. We could have let 
its global phase depend on time, for instance, and this would have given in Ea. (|26|) a contribution 
along \uj ) (\u> ) is defined in |Appendix A ). Similarly, a more general variational ansatz could 



also have been used in Appendix B[ The important point is that if the soliton's parameters are 



chosen within the same variational space, their time evolution is described - in the adiabatic and 
perturbative approach - by the same equations. Besides, the radiative term (pi having in all cases 
to be restricted to the phonon part of the spectrum, its time evolution is not (at least in the limit 
V 2 ^> e U; see below) affected by the specific choice of variational parameters used for describing 
the soliton. 

3.2 Radiated part 

The time evolution of the radiative component |<^>i) is obtained in a manner similar to what is done 
for the soliton's parameters. Projecting Eq. (|26[) onto the phonon eigenfunctions of TC by applying 
(S°|cr 3 yields 

iN« d t0 C« = Mq e a q C% + (E°\K) , (31) 
where a = ±. e° in Eq. (|3l"j) is the eigenvalue of TL associated with [see Eq. (jA2|) ] 

e° = q (^-sme + a^+l] , (32) 

and Mq is a normalization factor [see Eas. (|A4|) and (|A5[) ]. In deriving Ea. (j31|) . we have taken into 
account that the eigenfunctions depend on t only through the slow time ti (via sin#). The 
same holds for and Mg. Thus, writing 

C%(to,ti) = D%(to,ti) exp{-i e°to} , (33) 

one has, at the same order of approximation as Eq. pif) . 

1 .„ 



dt ^ = Tj^^ n) - (34) 

In integrating Eq. (|34[) we can choose between two equivalent strategies. The first (and difficult) one 
is to solve this equation taking into account that t\ = e to an d that 9 and z have the time dependence 
specified by Eas. (P3*)) .(|2^ |) and (|3U|). The second one is to integrate this equation considering to and 
t\ as independent variables. In this case, the t\ dependence of 9 and z will not matter and the to 
dependence of z will be specified by Ea. (|24l) . According to this second method one obtains 



Wo,<i) = f t0 dt' e i£ ^ (E°\K) + D-(t 1 ) , (35) 
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where Dg is an unknown function of ti [verifying D^(ti — * — oo) = 0] which could be determined 
by pushing the perturbative expansion to next order in e. In the following we will simply neglect 
this term. This is legitimate in the limit where all the ti-dependent terms are nearly constant-i.e., 
to the limit where the parameters of the soliton are very weakly affected by the obstacle. We will 
see in the next section that this limit is reached when eU <C V 2 . 



Of most interest to us is the total amount of radiation emitted by the soliton. For the deter- 
mination of this quantity we need the explicit expression of Eq.(J3SJ) at large times. In the limit to 
and t\ — > +00, Eq. (|35l) (without the D® term) reads explicitly 

1 



D«(+oo) 



dx dt' e 



1F a t' 



u,, 



'(x) TZ(x, x + z(t' , +00)) + v, 



(x,z(t' ,+oo)) , (36) 



where the functions u^(x) and v^{x) are the explicit components of |S°) defined in Appendix A 
[Eq. (|A3|) ]. Note that the ii-dependent parameters in Ea. (|3*6^1 have been given their asymptotic 
value. In particular, $m9(t\ — > +00) = sin#o = V, and according to Eq. (|24|) one has here 
z(t' , +00) = Vt' + z(+oo). The integration along t' Q in Eq. ()36j) can be computed easily using the 
expression (|2T|) for 1Z and Eq.® for Sf, leading to 



D«(+oo) 



-ie«2(+oo)/V 



dx e - i£ % x/v d x § 



+ <*(x) $(:r 



IV 

A long but straightforward computation gives the final result 

1 q U*(e^/V)e- l£ ^ ( - +Q °y v 



D"(+oo) 



16 ^1 + g 2/4 



sinh 



(37) 



(38) 



2 y vi - v 2 

In this formula the term z(+oo) can be obtained through the numerical determination of z(t). 
We indicate in section 14.21 different approximation schemes allowing one to obtain an analytical 
evaluation of this term [Eq. (|47|> and below]. ^From expression (|38jl we see that the radiation 
contributes to ()2(J|) to the total wave function with a contribution of order (e U /V 2 )y/l — V 2 . 
According to the approximation scheme defined in the beginning of the present section we have 
V 2 3> e U . Since U and U are of same order of magnitude, the radiated part is, as expected, a 
small quantity. 



4 Analysis of the results 

In this section we analyse the solutions of Eqs ()23 |) .([29 j) . and (|30[). I|28 jl .l|33j l and which describe 
the dynamics of the system within our approach. The separation between the slow and fast times 
we used up to now in order to identify which time derivatives were negligible is no longer necessary, 
and we will henceforth only employ the actual time t. We will also drop the multiplicative factor e 
in front of the perturbing potential U(z) and of 4>i(x, t). In the two following subsections we study 
the evolution of the parameters of the soliton and in the two last ones we analyse the radiated part. 
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4.1 Effective potential approximation 



Since we now use the actual time t, instead of using Eas.(|23 j) .(|29 [) and (|3U|) . it is more appropriate 
to work with the equivalent equations (|B9IB10|) , In order to get insight into the details of the 
dynamics of the soliton, one should solve these equations numerically for a particular obstacle. 
This is done in section 14.21 where we study the behavior of a soliton incident on a delta scatterer. 
But before going to this point, it is interesting to study some limiting cases. In particular, the 
dynamics of the variationnal solution (|16|) can be more easily understood in the limit of a very 
dark soliton (almost back). To this end, let us multiply Eq. l)B9j) by z and add it to Eq. QBlOj) 
multiplied by 9. This gives 

4 9 sin 9 cos 2 9 = 2 Re j J dz R* [cf) sol ] (z d^ sol + 9 <%0 sol ) j • (39) 

In the limit of a weak potential, we have to keep in mind that R is a small quantity (of order 
of U). It is then legitimate at first order to replace on the RHS of (|39() z by sin# and to drop the 
term 9. One thus obtains 



3 2f) f + °°, dj 

- cos 9 dz j- ■ 

4 J-oo cosh [cos 9{z — z)\ 



cos 2 9 j dz ZJ n/ — . (40) 



A further simplification of the equations is obtained in the limit of very dark soliton, when 
-> 0. In this limit 9 ~ z and using expression (JHJ) for / we can put Eq. (|4U|) in the following form: 



^^ ff where U ciS {z) = — - f dz — — = ~ [ dz ^ - 

dz 2 J _ OQ cosh (z — z) 2 y_ 00 cosh (z - 



If we furthermore consider a potential U (z) which slowly depends on z (over a length scale much 
larger than unity [38]), then U(z) in the convolution of the RHS of Ea.(|41[) does not appreciably 
vary over the distance where the term cosh _2 (z — z) is noticeable. This yields 

u(z) r +co i 

U cB (z) ^^r 1 dz —j- = U{z) . (42) 

I j-oo cosh (z — z) 

Equations (j41|) and 1)42(1 show than in the appropriate limit (very dark soliton, weak and slowly 
varying potential) the soliton can be considered as an effective classical particle of mass 2 (i.e., 
twice the mass of a bare particle) of position z (the position of the center of the density trough) 
evolving in a potential U(z). If we relax the hypothesis of slowly varying potential, the soliton 
can still be considered as a particle of mass 2, but it now evolves in an effective potential U eB (z) 
defined in Eq.(|41|) as a convolution of the real potential U(z). The fact that the effective mass of 
the soliton is twice the one of a bare particle has already been obtained in Refs. [16, 17, 19,39]. 
Previous studies mainly focused on slowly varying external potentials and, as a result, the existence 
of an effective potential U cS - different from U - had not been noticed so far, except in Ref. [17] 
where this result has already been obtained in the special case of a 5 scatterer. In the following, we 
denote the approximation corresponding to Eq. (|41|) as the effective potential approximation: the 
soliton is considered as an effective classical particle of mass 2, position z, moving in the potential 
U eS (z). 
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4.2 Numerical check 



Let us now study in detail a particular example. We consider a soliton incident on a pointlike 
obstacle-i.e., a 5 scatterer characterized by U(x) = X5(x). In this case, the static background f(z) 
is given by Eq© and Eqs. ()B9|) and (jBlOj) read 



sgn(A) cos 



dz 

sinh(2z + 2a) ^ oo7h 4 X cosh 4 Y , 
cos 3 6> f +QO r , 9/ ., /tanhX tanhYA 



and 



, / X cosh" 2 X + tanh X Y cosh" 2 Y + tanh Y 
suit) — z = sgn(A)sm# / . — ^ l~ 



n 



sinh(22; + 2a) \ cosh 2 X cosh 2 Y 



sin cos f + °°, r „ 9 , /1-XtanhX l-YtanhY\ , . 

+ L dz [1 " ; (z)] ( «*fx + J ' (44) 

X and y in Eqs. (|43|) and (|44|) are notations for (z — z) cos and (z + z) cos ^ respectively, and 
the expressions of function / and of parameter a are given in Eq. Solving Eqs. (|43|) and l|44|) 
numerically, we obtain the time evolution of the parameters of the soliton. We plot in Figs. 2 
and 3 the behavior of z as a function of t for different initial velocities V. Figure 2 corresponds 
to a repulsive interaction with A = +1 and Fig. 3 to an attractive one with A = —1. The initial 
conditions for the numerical integration of Eqs. (|43jl and (|44[) are taken to be z(t = 0) = — 10 
and z(t = 0) = sin[0(t = 0)] = V. Several curves are drawn, corresponding to several values of 
V . In the repulsive case (Fig. 2), three initial velocities have been chosen: V = 0.9, 0.707, and 
0.4. The value V = 0.9 corresponds to a fast soliton which is weakly perturbed by the barrier, 
the value V = 0.4 corresponds to a reflected soliton, and the value V = 0.707 is just below the 
value V = y A/2 which, according to the effective potential approximation (|41|). is the separatrix 
between transmission and reflexion (corresponding to V 2 = max {U eS (z)} = A/2). In the attractive 
case (Fig. 3) the curves are drawn in the cases V = 0.707, 0.4 and 0.3. In both figures, the solid 
lines correspond to the exact numerical solution of Eas. (|43|) and (|44j) and the dashed lines to the 
result of the effective potential approximation. 

We first remark that the case of a 5 scatterer we consider here is the worst possible for the effec- 
tive potential approximation and that this approximation is certainly more at ease with smoother 
potentials. However, it is interesting to note that the effective potential approximation, which 
could be thought as oversimplified, is often very good. The worst agreement occurs in the case of 
repulsive obstacle, near the separatrix (which is estimated by the effective potential approximation 
to occur in the case of Fig. 2 at V = l/y/2). As we will see below (Fig. 4), the effective potential 
approximation does not exactly predict the location of this separatrix whereas, in this region, the 
trajectories are strongly affected by small changes of the initial velocity V. This is the reason for 
the bad agreement of the result of the approximate method with the ones given by the numerical 
integration of Eqs. (|43[) and (|44|) for V = 0.707. However, it is surprising to note that the effective 
potential approximation is generically valid, even in the case where the soliton is far from being 
very dark: even the limit V — > 1 is very accurately described by this approximation on Figs. 2 and 
3. 
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FIG. 2. z(t) for solitons of initial velocity 
V, incident on a repulsive obstacle U(x) = 
X5(x) with A = +1. The solid lines corre- 
spond to the numerical solution of Eqs. i|43)) 
and ifHj) and the dashed lines to the effective 
potential approximation 



FIG. 3.: Same as Fig. 2 for solitons inci- 
dent on an attractive obstacle U(x) = X5(x) 
with A = — 1. The dashed lines correspon- 
ding to the effective potential approximation 
are hardly distinguishable from the solid lines 
which correspond to the numerical solution of 
Eqs. USD and pi) . 



In order to investigate more precisely the limit of large initial velocities V and to assess the 
validity of the effective potential approximation, let us now establish the form of Eqs. ()B9|) and 
(Tmnjl in the case of a weakly perturbed soliton. From the effective potential approximation, one 
infers that the soliton is weakly perturbed by the obstacle when its initial energy is large compared 
to the external potential U cff - i.e., in the regime V 2 3> U (since U cff and U are typically of same 
order of magnitude). This is confirmed by the numerical results presented on Figs. 2 and 3: the 
trajectory of the soliton is less modified for large V. In the extreme limit V 2 ^> U one may write 
e(t) = 6»o + 0(i) and z = Vt + A(t), with 6 < 9 and A < V. A has the meaning of a shift 
in position: it is the difference between the position of the center of the soliton in presence of the 
obstacle with the value it would have in absence of the obstacle. The perturbative versions of Eqs. 
(fB9l) and (fBlOl) read 

/+oo 
dz R*[$(z-Vt,do)]<& z (z-Vt,d ) (45) 
-oo 

and 

/+oo 
dz R* - Vt, 9 )} <S> e {z - Vt, 6 ) . (46) 
-oo 

^From these equations it is a simple matter to compute analytically the asymptotic expressions 
of the soliton parameter. One obtains - as expected - 0(+oo) = 0, and the asymptotic shift in 
position is 

A(+oo) = -U(0) , (47) 



13 



where U(0) = J^dxU(x). Equation (ji7|) for A(+oo) is an approximation (valid in the regime 
V 2 » U) of the exact result 

A(+oo) = lim {z(t) - Vt} . (48) 

Comparing definitions l|24jl and (|48[) we see that, since sin9(ti — * +00) = V, one has A(+oo) = 
z(+oo). In the case of a 5 scatterer, the exact value Q48|) was computed through numerical solution 
of Eqs. and (|44|). The result is displayed in Fig. 4 (thick solid curves) and compared with 

the approximate expression (|47j) (thin solid curves) and with the result of the effective potential 
approximation (dashed curves). 




FIG. 4. A(+oo) as a function of the initial 
velocity V of a soliton incident on a 5 peak 
U(x) = X5(x). The upper curves correspond 
to the case A = —0.5, the lower ones to the 
case A = 0.5. The thick solid lines are the 
exact result i)48J) obtained form the numeri- 
cal integration of Eqs. and ifHj) . The 
dashed curves are the result P^j) of the ef- 
fective potential approximation and the thin 
solid curves are the approximate result ipIT)) , 



In the case of the effective potential approximation, the value of the shift A(+oo) can be 
computed either via the numerical solution of the equation of motion (|41j) or via the formula 



A(+oo) 



dx 



y/1 ~ U ctt {x)/V 2 



(49) 



^From this expression, one sees that in the limit V 2 3> U cH ~ U, the effective potential approx- 
imation yields a result A(+oo) ~ — ^U(0)/V 2 . Hence, in this limit, the shift computed via the 
effective potential approximation at V = 1 is correct [since it agrees with the result (|47(l at V = 1]. 
This is surprising, because the effective potential approximation is expected to be accurate only for 
very dark solitons. However, one can also notice that detailed agreement with the exact result (|48|) 
is missed since, in the limit V 2 3> U, the asymptotic evaluation (JITj) of (Jl%|) does not exactly match 
the one of (|49|) . Yet one sees from Fig. 4 that the shift computed via the effective potential approx- 
imation is in surprisingly good agreement with the exact value, even for fast solitons. In particular, 
in the case of an attractive potential, the exact evaluation of A(+oo) and its approximation (|49[) 
are hardly distinguishable. 



4.3 Backward- and forward- emitted wave packets 

At this point it is interesting to study in more detail the structure of the phonon part of the 
wave function-i.e., of 4>i(x,t). ^From Eqs. and ()A3|) one can separate 4>\ into two parts: 
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+ <p 1 with 

tf(x,t)= ldqC«(t)u a Jx) . (-,()) 



iFvom the explicit expressions Q32JI . l|3^|) and ()A3j) one sees that <fi^ (0^) describes waves propa- 
gating toward the positive (negative) x. 

We are interested in studying the outcome of the collision-i.e., in obtaining an analytical eval- 
uation of Ea ()50|) when t — > +00. To this end, one uses the fact that at large time one has 
Cg(t) oc exp{— is® [t + z(+oo)/V]}. Hence, instead of working with the variable t, it is convenient 
here to define r = t + z(+oo)/V and to write Eq. (J5(J|) in the form 

<ft(x, t) = [dq G a (q, x) exp{i[g(x + Vr)-ar F(q)}} , (51) 

Jr 

where F(q) = q (q 2 /4 + I) 1 / 2 and G a {q,x) = [q/2 + e*/q + ix(x)] 2 D% (+00) exp{ie£z(+oo)/F}. In 
the appropriate limit (to be defined soon), one can evaluate this expression through a saddle phase 
estimate. In this limit, the rapidly oscillating phase in Ea. (|51|1 is stationary at point ±q a which 
are solutions of x + V r = a r F'(q). One has 

1 r„o „ rrTT, rl x 



2 



X 2 - 4 + a X Vx 2 + 8 



with X = V + - . (52) 



T 



One can easily verify that q a goes to zero when = a and that q\ is positive only if a(V+-) > 1. 
From this, one sees that the saddle phase estimate of Eas.([5U|) and (|51[) is accurate when the two 
saddles are well separated-i.e., in the regime x 3> (1 — V) r for a = + and x <C —(1 + V) r for 
a = — . If this condition is fulfilled, one obtains 



</,'(•'•■/> - (;"|>/ n ./) \l \Jj ^lq<,(x+Vr)-arF( ga )~an/4] 

r \Qa)\ 



+ G a {-q a , X ) x \ J* e -i[l^+Vr)-arF( qa )-an/A] _ ^ 
l-r [Qa)\ 

The exact expression computed from Ea. (|5Uj) is compared in Fig. 5 with the saddle phase estimate 
H53JI . The curves are drawn at r = 60 [40] for a soliton with incident velocity V = 0.5. The obstacle 
is here taken to be a delta scatterer \5(x). <f>\ being proportional to A [through the expression 
(|38p of Dg(+oo)] we represent in Fig. 5 the value of (pi(x,t)/X (actually its real part) which do no 
depend on A. 



One sees in Figure 5 that the semiclassical approximation (|53|) is excellent in all its expected 
domain of validity and diverges at x = (1 — V) r = 30 (for a = +) and x = —(1 + V) r = —90 
(for a = — ) [41]. Hence, these points can be considered as representative of the region where the 
contribution of 4>f and 4>± to the total wave function is more important. Roughly speaking, the 
present approach indicates that, long after the collision, (f>i(x, t) is maximum around x = {a — V)r. 
We recall that when using x (instead of z) as position coordinate, the soliton is, at all times, located 
around x = 0. Hence, going back to the z coordinate, we have a clear picture of the process at large 
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times: the soliton propagates at velocity V (the same as its initial velocity) after having emitted 
phonons which form two wave packets, one propagating in the forward direction with group velocity 
1 (i.e., the sound velocity) and the other one propagating backwards with group velocity —1. The 
same conclusion seems to be reached in the numerical simulations of Parker et al. [21,22]. 



Fig. 5. Re{(pi(x,t)} as a function of x 
for t = 60 for a soliton of initial velo- 
city V = 0.5 incident on a 5 scatterer 
X5(x). The thick line represents the re- 
sult i|50j) and the thin line its semiclas- 
sical approximation $5<1\) . For legibility 
we have separated the region where cft^ is 
nonzero (around x = —90 ) from the one 
where <pf is non zero (around x = 30). 
Note that in the expected domain of va- 
lidity ofEq. (x > 30 and x <C -90) 
one can hardly distinguish the thick line 
from its semiclassical approximation. 



4.4 Radiated energy 

A quantity of importance for characterizing the system is the total energy radiated by the soliton. 
Equation© for the field ip which, in the present work, is of the form (f(z, t) exp(— i t) [cf. Eq. (|10|)]. 
conserves the energy £ defined as: 

£[<p] = jdz {- \cp z \ 2 + \ (M 2 - l) 2 + U(z) M 2 } . (54) 

In order to have an expression of the energy in terms of the field <p which, when (f> = matches 
the usual expression (jB3|) of the energy of the soliton, we rather work with the quantity E[(p} = 
£[f(ft] — £[/]■ E[cp] is of course a conserved quantity, and we are interested in its expression far 
before (t — ► — oo) and far after (t — > +oo) the collision with the obstacle. We note here that f(z) — 1 
and U(z) are non zero only when z is close to to origin, whereas, in the same region, <p(z, t) — 1 is 
zero when t — * ±oo. After a change of variable from z to x = z — z(t), the previous remark allows 
one to obtain the simplified expression for E (only valid when t — ► ±oo): 

E[4>] = \ fdx{\4> x \ 2 + (\4>?-l) 2 } • (55) 

Using the decomposition (|15j) . keeping the lowest orders in 5<p, and taking into account the fact 
that, when t — > ±oo, |<1>| 2 — 1 is zero in the regions where Sep is noticeable, one obtains 

E[4>] = \ [dx {\$ x \ 2 + m 2 - l) 2 } + ~ fdx {|5^| 2 + (<r<50 + <5<f $) 2 ) + 0(S<p 3 ) . (56) 
2 JR 2 J R L ) 
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The first integral on the RHS of Eq. ()56|) corresponds to the soliton's energy and is equal to 3 cos 3 0. 
The second integral on the RHS of Eq. ()56|) corresponds to the energy of the radiated part and is 
denoted by E xad in the following. 

We are now facing a difficulty: we performed a computation at order e and at this order 
we have #(+00) = #(—00) since the equations for the parameters of the soliton are the same as 
the one obtained in the adiabatic approximation (see the discussion at the end of Appendix B ). 



Accordingly, E rad in Eq. (|5Hj) being of order e 2 should be neglected. Hence, at order e nothing has 
occurred for the soliton's energy: this quantity is not modified by the collision with the obstacle 
and the radiated energy should be neglected. Thus, it seems that our first-order approach is unable 
to predict the amount of energy lost by the soliton during the collision with the obstacle. 

However, as already remarked in the study of the scattering of bright solitons [42], one can cir- 
cumvent this difficulty and extract some second-order information from our results. The procedure 
is the following: when pushing the computations at order e 2 , the 0(e 2 ) estimate of E vad is still 
given by the second term on the RHS of (|56|) with 5(j> = e<pi, which we know from our first- order 
approach. At second order, since -E rad is non zero, the soliton's energy has been modified by the 
collision and energy conservation now reads 

I cos 3 [0(-oc)] = E = I cos 3 [0(+oc)] + £ rad . (57) 

Equation (|57jl allows us to determine the change in the soliton's parameter 0. Writing 0(— 00) = 9q 
(with sin = V) and 0(+oo) = 9q + 59 one obtains 

69 = -a — fr d - a • ( 58 ) 

4 cos 2 O sin O 

From Eq. (|58jl one can also determine the velocity at t — * +00 which is equal to sin[0(+oo)] = 
V + 59 cos 0o- Thus, we can determine how the collision has affected the soliton's shape and 
velocity by computing E iad (replacing 5(f) by 4>\). This will be done in the remaining of this section. 

On the basis of the analysis in terms of forward- and backward-emitted wave packet made in 
Section 1431 one can separate E iad into two parts, which we denote E~ d and E^ d , the first one 
corresponding to energy radiated backwards and the second one to forward-radiated energy, with 



i Jdx {|5^| 2 + ($*5(f> a + 5<f> a *$) 2 } . (59) 
A long computation which is summarized in Appendix C yields the result 



E a , = lim - 



EZ d = 167T £°°dq |^(+oo)| 2 (e«) 2 + ■ (60) 



When Dq(+oo) is given by Ea. ()38|) . one obtains 



7T 



„ , i 2 (e?) 2 \U(£?/V)\ 2 

e ™* = Wvg d i 7 rr= 4m-- ( 61 ) 
16 v*h sinh2 f^7TTg74 

2V VI - V 2 
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The behavior at low and high velocity of E® A defined in Eq. ()61|) is the following 

H-oo a iff, \i2 

E -wl whe " (62) 

and 

±(l-^ )5 /2|£ (0) |2 i ^^Ifl-ya^l^O)! 3 , when 7-1. (63) 

One sees from Eq. (|62|) that our approach predicts an unphysical divergence of the radiated energy 
at low incident soliton velocity. On the contrary, numerical computations indicate that a soliton 
with very low velocity does not radiate [21,22]. However, one must bear in mind that H61|) is the 
result of a first-order expansion only valid in the limit V 2 3> U and is unable to tackle the regime of 
very low incident velocities. More interestingly, in the high-velocity regime - where the first-order 
perturbation theory is valid - we see from Eq. (|63j) that the leading-order estimate of the total 
amount of radiation (forward or backward emitted) vanishes. 

In order to fix the ideas, we plot in Fig. 6 the value of E® ad as a function of the initial soliton 
velocity V. The obstacle is here taken to be a delta scatterer U(z) = XS(z). In this case U(q) = A. 



FIG 6. Energy Ef ad radiated in the for- 
ward (a = +) and backward (a = — ) 
directions by a soliton of initial velo- 
city V incident on a 5 scatterer. The 
solid lines represent the result If61)) and 
the dashed line the approximation 
which reads here E^ d ~ A 2 / (15 V). 
The inset displays a blowup of the fi- 
gure at high velocity In the inset, the 
dashed curves are the asymptotic re- 
sults JS3J). 

" 0.2 0.4 y 0.6 0.8 1 

Figure 6 shows that most of the energy is radiated backwards (this was already implicit in Fig. 
5) and confirms that, at leading order in U/V 2 , a soliton does not radiate in the limit V — * 1. 
Besides, not only the absolute value of E ia _ d goes to zero, but also the relative amount of energy 
radiated E rad /E vanishes [as (1 — V 2 )]. Very similar results are obtained for an obstacle interacting 
with the beam through a finite-range potential (for instance, a Gaussian). This absence of radiation 
of a fast soliton can be explained intuitively as follows: whatever the sign of the potential describing 
the obstacle, the soliton loses energy under the form of radiated phonons. Accordingly it gets less 
dark [56 > in Eq. (|58|) ] and is accelerated. This increased velocity after a loss of energy is 
a typical feature of dark solitons which are sometimes referred to as effective particles having a 
negative kinetic mass which decreases with increasing energy [43]. However, our results show that, 
since the soliton velocity cannot exceed the speed of sound, a soliton whose velocity is close to this 
upper limit cannot be further accelerated and the radiative process is suppressed. 
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5 Conclusion 



In this paper we have presented a study of the dynamics of a dark soliton experiencing a collision 
with a finite-size potential in a quasi-lD condensate. We determined the evolution of the soliton's 
parameters and also included radiative effects within secular perturbation theory. 

A first output of the present work is what we called the " effective potential theory" : in many 
instances the soliton can be described as an effective classical particle of mass 2 (twice the mass of 
a bare particle) evolving in an effective potential U eS [defined in Eq. (|41|)]. This approximation is 
rigorously valid in the case of a slow soliton incident on a weak potential, but its actual regime of 
validity appears to be quite broad. 

The effective potential theory is an approximation where - as in all adiabatic approaches - 
radiative effects are neglected. Perturbation theory allows one to get a deeper insight into the 
collisional process and to determine the amount of radiated energy at leading order in U/V 2 . We 
show that the radiated waves form two counterpropagating phonon wave packets, and we predict 
that the radiative process is suppressed in the limit of a soliton moving with a velocity close to the 
velocity of sound. This result should be checked numerically; work in this direction is in progress. 

Whereas adiabatic theory predicts that the soliton's shape and velocity are the same far before 
and far after the collision with the obstacle, it is an important feature of the perturbative approach 
of being able to determine finite asymptotic modifications of the soliton's parameters due to the 
collision. We computed [in Eq. (|58|) ] the modification of the soliton's parameters at leading order in 
U/V 2 . The qualitative picture of the collisional process drawn from our approach is the following: 
the soliton radiates energy, gets less dark, and is accelerated. Since the velocity of a dark soliton 
cannot exceed the velocity of sound in the system, it is natural that this velocity appears as a 
threshold for emission of radiations. Roughly speaking, a soliton with a velocity close to the 
velocity of sound cannot radiate [as seen from Eqs. (|6.3|) ] since its velocity cannot further increase. 
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Appendix A 

In this appendix we present the eigenvectors and eigenvalues of the Hamiltonian TC defined in 
Eq. (|27j) . H is not diagonalizable because, as we will see below, its null space and the one of H 2 are 
not identical. If we denote by its "generalized null space" [37] the union of these two null spaces, 
one can easily verify that it is spanned by the four vectors \oj ), \co e ), |^ ) ; an d |f2 e ) defined as 




= x + i sm $ 
— <J>* = — x + i sin# 



) 





(Al) 
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where the function x{ x i@) is defined in Ea. (|14|) , The kets defined in Ea, HAl|) verify Ti.\u) ) = 
H\u) e ) = and H 2 \VL ) =H 2 \Vt e ) = 0, with H\VL ) = 2 cos 2 9 \u) ) and H\tt e ) = cos 2 9 \u e ). One sees 
from Ea. (|Al|) that \uj e ) and |O e ) are, respectively, linked to variations of the center of the soliton 
and of the parameter 9 (i.e., to the phase change across the soliton): this is the reason why the 
terms in 9 tl and z tl in Eq. (|25|) can be rewritten in Eq. (|26|) by means of \ui e ) and |f2 e ). One can 
similarly show that \u ) is linked to modulations of the global phase of the soliton and that \Q D ) is 
linked to variations of the background density at infinity. 

The remainder of the spectrum of Ti is what we call the "phonon spectrum." It has two branches 
which we denote "+" and "— ." The corresponding eigenvectors and eigenvalues are denoted ) 
and with 

W|3±) = e±|3±). (A2) 

The explicit expression of the eigenvalues is given in the main text [Eq. 1)32(1 ], It can be simply 
obtained by considering the form of Eq. (|A2JI when x — > ±oo. In this limit, <I> goes to a constant, 
and looking for the eigenvectors under the form of plane waves, exjp{iqx} (JJq, V q ) T (where 
and V± are constants), yields the result (|32|). This is the reason why we denote these excitations as 
phonons. A better denomination should be "Bogoliubov excitations" because, far from the soliton, 
their form and dispersion relation correspond indeed to the elementary excitations of a constant 
background moving with velocity —V. 

The exact expression (valid for all x € M) of the eigenvectors is given by the squared Jost 
solutions of the inverse problem [23]. They read |H^) = (uf(x), v^(x)) T with 



u^{x) = exp{iqx} ( ^ + — + i% 



"q 

-± 
X 
2 q 



Vg(x) = exp{iqx} ^| - + ixj ■ (A3) 

The natural inner product of two kets is { - 1 c"3 1 ■ ) 5 where 03 is the third Pauli matrice. The 
eigenvectors have the following normalization: 

<^|o- 3 |H«) = f dx [uf(x)u«(x)-vl m (x)v%(x)] =Ar?6 ajP 6(p-q) , (A4) 

with 

A^ = 16avr g y?7T(l) 2 . (A5) 
In the main text we also use the following orthogonality relations: 

{uo e \a^\uj e ) = (O e |c73|^e) = (u) e \az\Eq) = (J7 e |a3|3g) = , (u e \a 3 \0, e ) = -4cos 3 9 (A6) 

and 

{aj \a 3 \oj ) = (uj \a 3 \uj e ) = (u \a 3 \Eg) = , (u \a s \£l e ) = 2 i sin 9 cos 9 . (A7) 
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Appendix B 



In this appendix we briefly present the Lagrangian approach for dark soliton of Kivshar and 
Krolikowski [14] and derive the Lagrange equations ()B9|) and (|B10|) , 

In absence of the perturbation R[(ft], Eq. (|11|) can be derived from the following Lagrangian 
density: 

C[cftA*\ = \ (<t>*<f>t - M*)(l - pa) - \\**? ~ \M? ~ !) 2 • (Bl) 
Accordingly, the energy and momentum are defined by 

dC 



E 



+00 

dz 
-00 
+00 

dz 

-00 



dcftt 



c 



dC 



+ 



"ft 
dC 



dz i\4> z \ 2 +m 



dz (<^*-<^> z )(l 



-i> 
1 

W\ 



(B2) 



The Lagrangian density IjBlj) is not a priori the most natural one leading to Eq. 1)11(1 . but for 
the asymptotic boundary condition we are working with (|0| — > 1 when z — ► ±00), it yields a 
finite value of the energy and, besides, the energy and momentum are now, for a field of the form 
<j)(x — Vt) (in particular, in the case of a soliton), related by the relation 5E = V SP, indicating that 
the background contribution has been removed and allowing one to treat the soliton as a classical 
particle-like object [33,44]. For completeness, we note that, for a soliton, (j) is given by Ea. (|13|) and 
its energy and momentum defined in Eq. (jB2|) have the following expressions: 



E = - cos" 



P = 7 r-26»-sin(2i 



(B3) 



Following Kivshar and Krolikowski [14], one can obtain adiabatic equations of motion for the 
soliton's parameters in the following way. Let us consider a variational approximation of the type of 
Eq. Q16|): the field of the soliton is parametrized with time dependent quantities q% (t), q n (t) and 
has no other time dependence: (f> so i(z, t) = (p(z, qi(t),..., q n (t))- One first defines the Lagrangian for 
the qiS as being 



L(q 1: qi,...,q n ,q n ) = 
Then the quantities d qi L and d^L are computed via 



+00 
dzC[ 

00 



(B4) 



and 



/+00 
dz (dq^d^C + d qi 4> z <90 z £ + d qi cj) t d^C) + c.c. 
-00 

/+00 
dz d qi 4> t d<p t C + c.c. 
-00 



(B5) 



(B6) 



where c.c. stands for complex conjugate. Considering that (ft is solution of Eq. 1)11(1 (including the 
perturbative term R[<ft]), simple manipulations allow one to obtain Lagrange-like equations for the 
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d f +OQ 
9 gi L - ^(<% L ) = J dz [d^C - d z {d^ z C) - d^d^C^d^ + c.c. 

= 2 Re { f °°dz R*^} d qi <f> sol ) . (B7) 



— oo 



In the particular case where 4> sol (z,t) = $(2 — z(t),9(t)) one obtains 



L(9, 9,z,z)=t[ir-26- sin(2 9)] - ~ cos 3 9 , (B8) 



and the equations of motion ()B7|) read explicitly 



r"+00 

4 9 cos 2 9 = 2 Re I I dz R* [<f> aol ] dgfa \ (B9) 



-00 



4 cos 2 9 (sm9 - z) = 2 Re \ / cte i?*[0 sol ] 5^ sol ^ . (BIO) 



and 



We note here a general feature, always valid in the framework of the adiabatic approximation: with 
equation (J2J) conserving energy, one can show that the soliton's energy defined in Eq. ([B2|) has the 
same value far before and far after the collision with the obstacle (the demonstration is essentially 
the same as the one given in Section T4. 41 where, in addition, the consequences of soliton's radiation - 
neglected in the present adiabatic approximation - are taken into account). As a result, one obtains 
for the solutions of Eqs ()B9|) and (|B10|) that #(+00) = 0(— 00), and l(±oo) = sin #(±00) = V. Hence 
the soliton's shape and velocity may change during the collision, but they eventually regain their 
initial values. This is intimately connected to the neglecting of radiative effects in the adiabatic 
approximation . 



Appendix C 

In this appendix we briefly indicate how to obtain expression for the radiated energy starting 
from Ea.(|59|l. where 5cj) is given by <f>i-i.e., by (|50|). Instead of giving a detained explanation on 
how to treat all the terms in the integrand of (|59|) . for brievity we focus on one of the contributions 
to the expression 1)59(1 for E^ d : 

f + OO f + OO / 2 Q \ 

/ dxm 2 \5^\ 2 = / dx 1- C 2 ° S |ff(M)| 2 . (Cl) 

J— 00 J— 00 \ cosn [x cos u j j 

We recall that we are interested of the evaluation of this term at large times. Expressing §\ through 
Eq. (|50|) . one can show that the term in cosh -2 in the integrand on the RHS of Eq. ()Cl|) can be 
dropped because it gives a contribution which decreases algebraically at t — ► +00 (this can be 
checked by a stationary phase evaluation of the integrals over the momenta). It thus remains to 
evaluate 

/+00 r+00 r+00 f+OO 

dx\^\ 2 = dq dp dx {C«C«\u«u«* -v«*v;*) + C«C; a v«v; a ) . (C2) 
-00 J—oo J —00 J —00 
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In Eq.,^, _ _ „ „ p 

normalization (|A4|) , For the evaluation of the last part of the integrand on the RHS of Ea. (|C2[) . 
the explicit expressions (|A3|) of u" (x) and (x) are to be used. In the course of this computation, 
an argument of stationary phase shows that only the x-independent terms with p = q give a finite 
contribution at t — > +00. These terms will contribute as 2ir5{p — q) after the integration over x. 
Altogether one obtains the expression 



dx \4>1 



a\2 



t— >+oo J _ 



dq\C, 



a\2 
1 I 



Q q \ , 2 
- - — I + COS ' 
2 q 



(C3) 



Noting that N„ defined in Eas HA4j) and (|A5|) is an odd function (and thus does not contribute to 
the integral since |C"| 2 is even) and explicitly computing the other contributions, one obtains 



dx |0fOM)f 



vr / dq (q 2 + 2) 



2 ' o,/^ 2 ljD « (+00) | 2 



Jo \ q 

The others contributions to Eg. ((!)§]) can be computed similarly. One obtains 

.o\ 2 



dx \d x <t>?\ 2 ~ 8vr / dqq z (q 2 + 2) -2- |D„ Q (+oo) 
Jo V Q 



and 



+00 



dx 



16 vr 



+00 



-Q\ 2 



Gathering all these contributions yields the final result 



(C4) 



(C5) 



(C6) 
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